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摘要 :【 AMY) sie Bombyx mori MX} ABE Beauveria bassiana 侵 染 的 应 管 基因 ,以 进一步 研究 家 看 抵御 真菌 
侵 染 的 分 子 机 制 。[ 方 法 ] 采 用 新 一 代 Solexa 高 通 量 测序 技术 对 感染 日 僵 菌 及 未 感染 白 僵 菌 的 对 照 组 家 午 进 行 了 
测序 分 析 , 得 选 差异 表达 基因 ; 结合 生物 信息 学 工具 分 析 差 异 表达 基因 的 功能 注释 、 分 类 及 涉及 的 信号 通路 等 ; 
应 用 菊 光 定量 PCR 技术 验证 10 个 基因 的 差异 表达 。[【 结 果 】 通 过 测序 和 生物 信息 学 分 析 共 获得 377 个 差异 表达 基 
， 其 中 表达 上 调 基因 236 +, 下 调 基因 141 个 ; KEGG 通路 分 析 表 明 , 各 通路 中 既 有 表达 上 调 的 基因 , 也 有 下 调 
基因 ; 12 个 上 调 基 因 、26 个 下 调 基 因 参 与 3 个 显著 性 富 集 的 KECC 通路 , 即 核糖 体 、 氨 酰 tRNA 生物 合成 和 剪接 
体 通路 。 定 量 PCR 与 测序 结果 显示 , 溶菌 酶 、 热 激 和 蛋白 、 谷 胱 甘 肽 S$- 转 移 酶 、 肽 聚 糖 识别 蛋白 等 与 免疫 应 激 相关 
的 蛋白 基因 均 呈 现 表 达 上 调 。【 结 论 】 本 研究 筛选 获得 的 差异 表达 基因 , 特别 是 上 调 表达 的 基因 可 能 与 家 和 蛋 应 对 
日 僵 菌 侵 染 的 应 答 机 制 有 关 ,， EPR REMO AE AAA AA AE A SEHR $- 转 移 酶 、 肽 聚 糖 
识别 重 白 基 因 等 可 能 直接 参与 了 家 看 对 白 伪 菌 的 免疫 识别 和 防御 , 研究 结果 为 从 分 子 水 平 阐明 家 看 抵御 真菌 侵 染 
的 防御 机 制 和 日 僵 基 对 家 和 蛋 的 致 病 机 理 提 供 新 的 依据 。 
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Differentially expressed genes in the silkworm (Bombyx mori) against 


infection by Beauveria bassiana 

HOU Cheng-Xiang ^ , QIN Guang-Xing ^ , GENG Tao'”, GAO Kun'”, PAN Zhong-Hua , QIAN He- 
Ying ^, GUO XiJie"^" (1. Sericultural Research Institute, Jiangsu University of Science and 
Technology, Zhenjiang 212018, Jiangsu, China; 2. Sericultural Research Institute, Chinese Academy of 
Agricultural Sciences, Zhenjiang, Jiangsu 212018, China) 

Abstract: [ Aim] To screen the responsive genes of the silkworm (Bombyx mori) against infection by 
Beauveria bassiana in order to further elaborate its molecular mechanisms of defensing B. bassiana 
infection. [ Methods] High throughput Solexa sequencing technology was employed in the transcriptome 
analysis of silkworm larvae of the infected and uninfected ( control) groups. Differentially expressed genes 
(DEGs) and their functions, classifications and pathways were analyzed by bioinformatics methods. Ten 
DEGs were further identified by fluorescence quantitative PCR. [Results] With Solexa sequencing and 
bioinformatic analysis, 377 DEGs including 236 up-regulated genes and 141 down-regulated genes were 
identified. KEGG pathway analysis indicated that there were both up- and down-regulated DEGs in any 
pathway. Twelve up-regulated DEGs and 26 down-regulated ones were involved in three significantly 
enriched KEGG ( Kyoto Encyclopedia of Genes and Genomes) pathways, i. e., ribosome, aminoacyl- 
tRNA biosynthesis and spliceosome pathways. Both the quantitative PCR and Solexa sequencing results 
indicated that the immune response-related genes including lysozyme, heat shock protein, glutathione S- 
transferase and peptidoglycan recognition protein genes were all up-regulated. [ Conclusion] The 
differentially expressed genes identified in this study, especially those up-regulated ones, are possibly 
involved in the silkworm response against B. bassiana invasion. The genes related to immune response 
such as those encoding lysozyme, heat shock protein, glutathione S-transferase and peptidoglycan 
recognition protein may be directly involved in immune recognition and defense of the silkworm against 


the B. bassiana infection. The results provide a new basis and clues for clarifying the mechanisms of 
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interaction between the silkworm and entomopathogenic fungi at the molecular level. 


Key words: Bombyx mori; Beauveria bassiana; 


expressed gene 


KÆ Bombyx mori 既是 重要 的 经 济 昆虫 , 也 是 
鲜 翅 目 昆 虫 的 模式 生物 。 真 再 病 是 牌 业 生 产 中 的 四 
大 传染 病 之 一 , ABE Beaueria bassiana 既是 其 中 
危害 最 严重 的 一 种 真菌 , 也 是 农林 害虫 生物 防治 中 
RAR AA BAC. ARRAS ECA 
相互 作用 是 和 宿主 感染 病原 并 激发 免疫 防御 的 生物 学 
基础 , 也 是 和 宿主 应 党 基 因 表 达 水 平 变化 的 起 因 。 日 
伪 菌 与 马尾 松 毛 虫 Dendrolimus punctatus (Lu et al., 
2008) . EM XE KM Ostrinia nubilalis (WERE, 
2011) 及 桑 天 牛 Apriona germari ( EjA%, 2009) 等 
昆虫 间 的 互 作 有 一 定 的 研究 , (ASA Ae BE OTE 
不 多 , 根据 相关 的 研究 结 采 难 以 建立 家 看 抗 真菌 的 
免疫 应 答 机 制 。 因 此 , RA a S 2C E ELE 
的 研究 ,寻找 宿主 体内 与 白 僵 菌 感染 密切 相关 的 应 
ZEN, 探索 其 作用 机 理 、 曾 明 其 表达 机 制 , 可 为 
家 至 真菌 病 的 防治 以 及 真菌 在 生物 防治 中 的 应 用 英 
定 基 础 。 

HEAR) HAS (fungi imperfecti) ,其 发 育 
周期 分 为 分 生 孢 子 、 和 营养 菌 丝 和 气 生 亢 丝 3 个 阶 
段 。 在 适宜 的 温 湿 度 下 ,附着 在 昆虫 体 壁 上 的 分 生 
孢子 侵入 得 主体 胁 并 伸 长 转化 为 能 吸收 养分 、 生 长 
繁殖 、 随 体液 扩展 到 全 号 的 营养 菌 丝 , VAS ZZ YS BE 
形成 能 脱离 母 菌 丝 、 目 行 吸 收养 分 、 同 一 端 或 两 站 
延伸 为 新 彰 养 苗 丝 的 茸 生 孢子 。 箱 主 死亡 1~24d 
后 ,， 菌 丝 穿 出 体 壁 形成 气 生 瑚 丝 , ZZ IE fl 
FEE RAIA De CBT EA. 

本 研究 将 Solexa 高 通 量 测序 技术 运用 到 感染 日 
僵 亢 的 家 看 转录 组 分 析 中 , 依据 家 看 基因 组 数据 
库 , 从 功能 基因 组 角度 分 析 家 看 应 对 白 僵 菌 侵 染 的 
ERA, Manu ot FKP ESP RR RS ER 
HE FE A E Pa SS B) B5 RS E] UT BLE fr EH o 


1 材料 与 方法 


1.1 材料 

本 实验 所 用 家 看 品种 为 本 所 保存 的 大 造 , 在 正 
fe ZR ARTF POST, RH 80% ~ 85% ) 饲养 至 3 龄 
x ( SIRS WIKRE 3 龄 幼虫 ) 后 分 为 2 组， 
每 组 60 3-FH-T- FASES S HIEN B. bassiana AH 
本 实验 室 从 看 体 上 复壮 、 分离 纯 化 并 保存 。 


immune response; transcriptome; differentially 


1.2 AER RR 

FAK AA RHE E FE E We i 1B sl ot HE fül 
F, 并 稀释 至 4 x 10°fl-F/mL, 346—428 E53 EK 
浮 液 中 浸泡 10 s 作为 感染 组 , 另 一 组 作为 对 照 在 灭 
菌 水 中 处理 相同 时 间 ; 感染 组 和 对 照 组 的 家 和 看 幼虫 
在 眠 期 已 用 看 座 消毒 粉剂 进行 体 表 消毒 。 感 染 处 理 
后 , 两 组 材料 均 在 28C RH 95% 的 高 温 高 湿 环 境 下 
饲养 8 h, 再 移 至 正常 温 湿 度 环境 中 饲养 。 接 种 8 h 
后 , 收集 感染 组 和 对 照 组 幼虫 各 15 头 ,， 各 平均 分 3 
I 迅速 投入 到 液 所 中 速 诛 ，- 80Y 保存 。 其 余 
的 看 继续 饲养 (正常 饲养 环境 下 ), 观察 生长 发 育 及 
发 病情 况 。 
1.3 总 RNA 的 提取 、 纯 化 与 cDNA 的 合成 

取 5 KA He oy, WE Ja Fl Trizol 试剂 
(Invitrogen 公司 , 美国 ) 提取 总 RNA。 定 量 后 取 
6 ug 总 RNA, Hj Illumina Gene Expression Sample 
Prep Kit (Illumina 公司 , 美国 ) 中 的 Oligo( dT) REER 
纯化 mRNA; 定性 、 定 量 后 用 PrimeScript® 1st 
Strand cDNA Synthesis Kit (TaKaRa 公司 , HÆ), 
按照 使 用 说 明 书 操作 步骤 反 转 录 合 成 cDNA, KE 
成 的 cDNA 样品 送 往 深圳 华 大 基因 公司 构建 数字 表 
达 谱 文库 。 
1.4 数字 基因 表达 谱 (digital gene expression 
profiling, DGE) 文库 的 构建 

利用 Illumina Cluster Station 和 Illumina HiSeq”” 
2000 系统 对 cDNA 样品 进行 数学 基因 表达 谱 文 库 
的 构建 工作 。 实 验 利用 Illumina Hiseq ^ 2000 系统 、 
采用 边 合成 边 测序 法 进行 测序 。 最 终 每 个 通道 将 产 
生 数 百 万 条 原始 图 像 数据 , 这 些 数据 经 base calling 
转化 为 原始 序列 数据 (raw data) o 
1.5 原始 序列 数据 分 析 

去 除 原始 序列 数据 中 I RR, MAUER SZ 
的 成 分 , 得 到 洁净 标签 (clean tag). ERKAM 
FE (ftp ://silkdb. org) 中 mRNA 上 所 有 的 CATG 位 
Ho 生成 CATG + 17nt 碱 基 的 参考 标签 数据 库 。 所 
有 的 clean tag 与 参考 标签 数据 库 进行 比 对 (最 多 人 允 
许 一 个 碱 基 错 配 ), 对 其 中 只 比 对 到 一 个 基因 的 标 
签 进行 基因 注释 ; 统计 每 个 基因 对 应 的 原始 clean 
tag ZU, 然后 对 原始 clean tag 数 做 标准 化 处 理 , 获得 
标准 化 的 基因 表达 量 。 


1 期 修成 香 等 : ERR ER EL HEA KRAMA ERA E A HT 15 


基因 注释 与 其 涉及 到 的 途径 分 析 主 要 参考 家 看 
数据 库 (ftp://silkdb. org) 和 KEGG 数据 库 (http:// 
www. genome. jp/kegg/pathway. html); 基因 分 类 根 
据 基因 本 体 论 (gene ontology，G0O) 中 细胞 成 分 
(cellular component) 、 分 子 功能 (molecular function ) 
和 生物 过 程 (biological process) 3 个 亚 类 进行 分 析 。 


1.6 RATE PCR 验证 分 析 

利用 Primer Premier 5.0 软件 设计 了 10 XJK 
定量 PCR 引物 , 引物 序列 及 相关 信息 见 表 1。 挑 选 
的 基因 中 部 分 可 能 与 昆虫 免疫 应 答 相 关 ， 如 DG3, 
DG5, DG7, DG8, DG15 和 DC16 ， 部 分 基因 参与 了 
代谢 、 物 质 转运 等 生理 过 程 , 如 DG2, DG7, DG11 
和 DG14; 而 DG9 则 参与 了 应 激 反应 。 


表 1 实时 谈 光 定量 PCR 引物 
Table 1 Primers used for real-time quantitative PCR 


基因 产物 


Gene product 


引物 名 称 


Primer name 


引物 序列 (5 -3") 


Primer sequence 


GTCTGCTACACGCTAAACTGCTCGC 


DG2 Hypothetical protein KGM_10198 
TTCCTTTTGGGCTTGTCTGTCTCTG 
TCAGCGACCACGAAATCAAAGAGG 
DG3 Cytochrome P450 CYP4G48 
TACCCAACAAACAGAGGACAAAGC 
CCTACCACAAGGCTTTAGATTTC 
DG5 Glutathione S-transferase omega 1 
ATCAGCGTATCCAGGTTCACTTC 
TCAGTTCACCCTCGCTTTTGGC 
DG7 Serine protease precursor 
TTTTGTTGCTTGTTGGCTCCCG 
GCGTGTGTGTGTTTGTGGACTTGTG 
DG8 Putative lysozyme 
TCTTTTATGTCTTCGTTTCGGAGCG 
AGGCAAACACGAGGAGAAGA 
DG9 Heat shock protein 19.9 
TGGGTGCGATTACAGACAAC 
TCAGCAGGGTGTCGAGATG 
DG11 Importin-7 
CGAGGGGTGTTGTGTAGCA 
AAACAAACTGTGATGTCGAAGC 
DG14 Asparagine synthetase 
CAGAGGTCTAACGCCATAAGGA 
DCIS Peptidoglycan recognition protein S6 AATTCCTTAGGCTGGGGTGA 
precursor CCGTGGACAGAAGCTTTTTT 
AGGATGGGAAGCAGAAGAAG 
DG16 30kP protease A (43k peptide) precursor 
GCGAGTAGCAGTGCAATGAG 
AATGGCTCCGGTATGTGC 
B-actin 
TTGCTCTGTGCCTCGTCT 


按照 SYBR Premix Ex Taq” Kit (TaKaRa 公司 ， 
HÆ) 说 明 书 进行 定量 PCR 实验 。 实 验 以 接种 处 
理 后 的 总 RNA 反 转 录 的 cDNA 为 模板 , 20 uL 反应 
体系 ; 反应 程序 为 : 95°C30 s; 95%C5 s, 60%C31 s, 
45 个 循环 ; 3 KERo URE b-actin 基因 为 内 参 ， 
用 2 3% ( Livak and Schmittgen, 2001 ) 分 析 家 看 
感染 日 僵 菌 后 感染 、 对 照 组 中 基因 相对 转录 水 平 的 
变化 。 每 处 理 和 对 照 重复 3 次 。 
1.7 数据 统计 分 析 

定量 PCR 实验 后 , 首先 分 别 算 出 每 个 生物 学 


感染 实验 中 目的 基因 与 内 参 在 感染 组 与 对 照 组 中 的 
3 个 C, 值 的 平均 值 , 然后 计算 感染 组 、 对 照 组 中 目 
的 基因 与 内 参 C, 值 的 平均 值 的 差 ; 再 以 前 者 减 去 
后 者 得 到 - AAC, (A; 最 后 计算 2 ^77, 所 得 的 数值 
就 是 目的 基因 在 感染 组 与 对 照 组 中 的 相对 表达 量 。 


2 结果 与 分 析 


2.1 Aerie Ae FRA RRR 
在 适宜 的 温 湿度 下 ,附着 在 家 看 幼虫 体 壁 上 的 
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Hf EAE 6 ~ 8 h 内 萌发 并 穿 透 体 壁 , 侵入 
昆虫 体内 。 接 种 处 理 后 约 36 h, 3 龄 幼虫 体 表 出 现 
TERRA ADE EDADES 感染 幼 
虫 在 60 h 内 全 部 死亡 ,显微镜 下 观察 濒 死 一 的 体 
液 , EMERE DERA; 约 72 h 后 死亡 的 幼虫 
体 壁 被 日 色 菌 丝 履 盖 。 这 说 明 处 理 乍 死 于 日 僵 亩 感 
染 ; 此 时 对 照看 处 于 正常 的 生长 状态 中 , 镜 检 体 液 
TORE fa ZZ 
2.2 ”数字 基因 表达 谱 文库 

分 别 构建 了 接种 8 h 后 感染 组 和 对 照 组 的 数字 
基因 表达 谱 文 库 。 原 始 图 像 经 转化 后 , 感染 组 与 对 
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照 组 库 中 分 别 得 到 3 969 200 和 4 819 322 个 原始 数 
据 ; 去 除 其 中 的 3 接头 、 低 质量 及 含 杂 质 的 成 分 
后 , 分 别 得 到 3 794 277 和 4 588 715 个 clean tag; 
其 中 , 29% ~31% 的 clean tag, 18% ~ 20% 的 Distinct 
clean tag 能 够 匹配 到 家 看 基因 组 数据 库 的 基因 中 。 
最 终 , 两 个 库 中 分 别 得 到 5 444 和 5 348 个 基因 , 其 
中 , 在 两 个 库 中 均 表 达 的 基因 有 4 034 ^P, 占 库 中 
基因 总 量 的 73.7% ; 仅 在 感染 组 中 表达 的 有 760 
个 , 713.9% ; 只 在 对 照 组 中 表达 的 有 678 +, 占 
12.4% 。 对 两 个 库 中 的 基因 比 对 、 汇 总 后 得 到 
5 472 个 不 同 的 基因 ( 表 2, 图 1)。 


R2 ， 家蚕 幼虫 白 僵 菌 感染 组 和 对 照 组 库 中 标签 (tags ) 与 基因 (genes) 的 分 布 


Table 2 Distribution of tags and genes in the Beauveria bassiana-infected and uninfected 


(control) libraries of Bombyx mori larvae 


标签 与 基因 


Tags and genes 


Ly 
BB 3 969 200 
Raw data 
洁净 标签 
TED 3 794 277 
Clean tags 
与 基因 匹配 的 全 部 标签 
3 BU RR 1 197 135 
All tags mapping to genes 
与 基因 匹配 的 明确 标签 
3 BEES 1 076 803 
Unambiguous tags mapping to genes 
与 标签 匹配 的 全 部 基因 5715 
All tag-mapped genes 
与 标签 匹配 的 明确 基因 5 444 
Unambiguous tag-mapped genes 
与 基因 组 匹配 的 标签 
3 Bei 2 026 647 
Tags mapping to genome 

NE 
RAN 570 495 


感染 组 
Infected group 


标签 总 数 独特 标签 所 占 百 分 比 
Total number of tags Percentage of distinct tags 
对 照 组 感染 组 对 照 组 


Control group Infected group Control group 


4819 322 4.58 4.20 
4 588 715 1.88 1.64 
1 338 847 20.26 18. 87 
1 215 887 19.56 18.21 
5 615 
5 348 
2 610 102 53.60 52.95 
639 766 26.15 28.18 


Unknown tags 



















Ries LI 
EN INN 







Infected group 4 034 Control group 
760 (73.7%) 678 
(13.9%) (12.4%) 


图 1 KAJE fe a BRS RIO RE PAE 
百分比 的 分 布 图 
Fig. 1 Distribution of number and percentage of genes in 


the Beauveria bassiana-infected and uninfected 


(control) libraries of Bombyx mori larvae 


2.3 GO 功能 注释 

根据 GO 分 类 中 分 子 功能 亚 类 分 析 本 实验 汇总 
后 所 获得 的 5 472 个 基因 , 发 现 有 885 个 基因 涉及 
了 90 种 分 子 功能 ; 由 于 一 个 基因 可 能 具有 多 个 功 
能 域 或 活性 ,因此 本 研究 以 每 个 功能 注释 下 的 二 级 
子 目 录 进 行 统计 分 析 。 统 计 结果 显示 : 排列 在 前 三 
的 分 别 是 具有 结合 (binding)、 水 解 酶 活性 
(hydrolase activity ) 和 结构 分 子 活 性 (structural 
molecule activity) 功能 的 基因 (图 2)。 
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X4 白 僵 菌 感染 8 h 后 家 等 体内 参与 KEGG 通路 的 基因 


Table 4 Genes involved in KEGG Pathways of Bombyx mori at 8 h post-inoculation of Beauveria bassiana 


KEGG 通路 016 
KEGG pathway Q-value 
FORE 0. 0000 
Ribosome 
tRNA 生物 合 
ARE Map 0. 0001 
Aminoacyl-tRNA biosynthesis 
Bt 
NER 0. 0003 
Spliceosome 
RARI 
色 氨 酸 代谢 0.2278 
Tryptophan metabolism 
EHE 0.2420 
Proteasome 
y - y = iA 
代谢 途径 0.4197 
Metabolic pathways 
BEE WOT 
ERRATA 0.4127 
Lipoic acid metabolism 
LJ Ab BE RT AN 
DURAESEEGSUN 0. 4127 
Antigen processing and presentation 
WALZ S S 
化 合 物 的 代谢 0 5233 
Selenocompound metabolism 
HAB; 
REUS 0. 5233 
Amoebiasis 
合计 
Total 


2.6 ”筛选 到 的 免疫 相关 基因 

病原 侵 当 昆虫 后 , 宿主 通过 模式 识别 受 体 识别 
病原 相关 分 子 模式 从 而 激活 昆虫 的 天 然 免 疫 , 而 胞 
外 级 联 的 重 日 酶 参与 调控 免疫 信号 的 扩大 或 缩小 ， 
再 由 具有 信号 转 导 功能 的 蛋 日 识别 、 转 化 免疫 信 
号 ,激发 宿主 分 泌 免 疫 因子 攻击 外 源 物 。 在 这 过 程 
F, 涉及 到 一 些 免 疫 相关 基因 。 本 实验 构建 的 文库 
中 有 一 些 和 表达 水 平 发 生变 化 的 免疫 相关 基因 C 
5) ,分 别 如 下 : 

模式 识别 受 体 : 本 实验 中 检测 到 多 个 能 编码 具 
有 识别 或 与 模式 识别 受 体 功 能 的 基因 ， 如 能 编码 
PGRP 56 precursor, B GRP 3 precursor, B GRP 4 
precursor, C 类 凝集 素 11 Hi [x (C-type lectin 11 
precursor). 、 分 子 量 为 50 kD 的 凝集 (50 kD 
lectin) ~ X i A (hemocytin), FA Wi HRW HA 3 
( thioester-containing protein, TEP). B 2S TÉ 3H KS 
(KA ( scavenger receptor class B) 等 蛋白 的 基因 。 

与 信号 调节 有 关 的 基因 : 本 实验 中 检测 到 发 卡 


P (H 基因 个 数 Number of genes 
P-value 上 调 Up-regulation 下调 Down-regulation 
0. 0000 2 12 
0. 0066 4 5 
0.0159 6 9 
0. 0053 3 3 
0.0071 1 4 
0.0174 
0.0177 1 1 
0.0193 4 0 
0. 0291 1 1 
0. 0306 7 1 

29 36 


ZETA Ia ZZ ARE A E (clip domain serine proteases, 
CLIPs). ARE A 8305) 3 Hi AK ( Serpin 3 前 
体 )、 丝 氨 酸 蛋 日 酶 抑制 剂 5 前 体 (Serpin 5 前 体 )、 
Serpin 13 前 体 、Serpin 14 前 体 、Serpin 21. Serpin 28 
等 与 免疫 信号 调 市 有 关 的 和 蛋 日 基 因 。 

与 信号 转 导 有 关 的 基因 : 本 实验 中 检测 到 Toll- 
like receptor 18-Wheeler, Toll 类 似 和 蛋白 (Toll-like 
protein). BE tF HU ^r 45 E A ( myeloid differential 
primary response protein, MyD88) ) Cactus, YM] 
Ed B WA (dorsal isoform B) €517 SE SKN 
H&R. 

免疫 因子 基因 : ANSE i A EE E] ENTE 
(putative lysozyme), X W HE Bj {X ( lysozyme 
precursor)、 酚 氧化 酶 激酶 前 体 (prophenoloxidase 
activating enzyme precursor), 2 fp WE PR SA 
(hemolin) , K## B HINA ( cecropin B precursor), 
抗菌 肽 (moricin I, lebocin 3 前 体 ) 等 免疫 因子 基因 。 
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RS 感染 白 僵 菌 的 家 等 数字 基因 表达 谱 文 库 中 表达 上 调 的 部 分 免疫 相关 基因 


Table 5 Up-regulated genes related with immunity in the DGE library of Bombyx mori infected with Beauveria bassiana 








xpi 感染 组 表达 量 对 照 组 表达 量 表达 量 比值 (感染 组 /对 照 组 ) 
C Expression level of Expression level of Ratio of expression level 
enes 
the infected group the control group ( infected/control ) 
糖 识别 蛋白 So 前 体 基因 
BIORBEUUILALET S6 前 体 54. 56 21.14 2.60 
Peptidoglycan recognition protein S6 precursor gene 
1,3 A REBELDE H 4 前 体 基 因 
B1, 5 RAEE 4 VE 697.10 215.53 3.20 
B-1, 3-Glucan recognition protein 4 precursor gene 
Ff 0. 53 0.01 52.10 
Hemocytin gene 
50 kD RIK] 
SEIS 145. 75 36. 18 4. 00 
50 kD Lectin gene 
B 类 清道 夫 受 体 3 基因 
HIERE 1.05 0.65 1.60 
Scavenger receptor class B member 3 gene 
££ 21148 Al J78] 3 前 体 基 因 
BURRIS BETHEL 3 HUE 19.77 2. 40 8.00 
Serine protease inhibitor 3 precursor gene 
丝氨酸 139] 21 基因 
AERIS RR 13. 70 3.92 3.50 
Serine protease inhibitor 21 gene 
22 ANG 13 基因 
AREAN 3.16 0.65 4.90 
Serine protease 13 gene 
HAN 前 体 基因 
SUP TA 5.80 2.62 2.10 
Serine protease precursor gene 
Cactus 基因 
aoas 6.85 4.14 1.60 
Cactus gene 
定 的 Toll 样 受 体 基 因 
BREI Toll SER 0.46 0.01 45.30 
Putative Toll-like receptor gene 
Toll 样 受 体 18-Wheel 
OIL PESEE IB Wheeler AEP 4.48 2.18 2.10 
Toll-like receptor 18-Wheeler gene 
EROR IEEE 
FOE ea 45. 59 14. 38 3.20 
Putative lysozyme gene 
溶菌 酶 前 体 基 因 
BEIM 65. 89 30. 29 2.10 
Lysozyme precursor gene 
Lea HK I 基因 
PLE 9.22 2.18 4.30 
Moricin | gene 
ER B HTA 48. 23 33. 78 1.40 
Cecropin-B precursor gene 
2.7 ”荧光 定量 PCR 结果 分 析 个 基因 的 相对 表达 趋势 与 数字 基因 表达 谱 的 实验 结 


挑选 了 10 个 差异 表达 的 基因 , 这 些 基 因 包 括 RM, BER T 30kP protease A (43k peptide) 
昆虫 拟 免疫 相关 基因 和 部 分 参与 代谢 、 物 质 转 运 等 precursor 基因 表达 下 调 外 ,其 余 基因 表达 均 上 调 
生理 过 程 的 基因 。 根 据 其 序列 设计 定量 PCR 引物 (图 5); 男 外 , 编码 PGRP S6 precursor 的 基因 差异 
进行 灾 光 定量 PCR 分 析 ( 表 1) 。 绪 条 表明 : 这 10 REAR G). 
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3 讨论 


FEAR ER LAA P, AEE BO ESE. MR 
急 的 一 种 病原 。 在 适宜 的 温 湿度 下 ,附着 在 看 体 壁 
上 的 孢子 可 膨大 、 穿 透 体 壁 侵入 宿主 体内 ; 随后 利 
用 宿主 养分 生长 、 繁 殖 , 最终 导致 宿主 死亡 。 在 整 
个 感染 过 程 中 , 与 孢子 萌发 及 其 在 和 宿主 体内 生长 繁 
殖 的 时 间 相 比较 , 穿 透 窒 主体 壁 所 需要 的 时 间 较 
短 , 但 其 感染 及 宿主 的 抗 性 机 制 依然 所 知 其 少 。 

本 人 研究 利用 转录 组 测序 技术 研究 家 看 锌 日 僵 落 
侵 染 初期 基因 的 差异 表达 ,发 现 一 些 具 有 识别 、 结 
合 、 信 号 调 市 与 转 导 、 人 免疫 因子 功能 的 免疫 相关 基 
因 表达 上 调 。 在 具有 识别 功能 的 蛋白 中 ，C-type 
lectin 11 是 家 看 能 特异 识别 真菌 的 企 日 ; hemocytin 
是 家 簿 特有 的 一 种 C-type lectin EH, 目前 仅 在 家 
等 和 采 蝇 等 昆虫 中 发 现 。 在 免疫 因子 中 ，hemolin 
和 moricin 是 家 看 特有 的 人 免疫 相关 和 集 日 , 日 前 在 果 
晶 、 按 蚊 等 昆虫 中 尚未 发 现 (Tanaka et al., 2008) 。 
与 抑制 性 消减 杂交 与 世族 等 方法 相 比 较 ，Solexa W 
序 方法 快捷 、 精 确 , 且 获 得 的 差异 表达 基因 不 仅 数 
ES, AY RMP) (Xiao et al., 2010; Hou et al., 
2011; Wu et al., 2011) , 

本 研究 的 KEGG 通路 分 析 表 明 , 共有 38 个 差 
异 表达 基因 参与 了 3 个 显著 性 寅 集 的 途径 ,分别 是 
核糖 体 、 氮 酰 -tRNA 生物 合成 及 剪接 体 通 路 。 日 僵 
P181 AC, 其 体内 的 负 丝 既 要 “躲避 ”和 宿主 
人 免疫 系统 的 识别 ,又 要 繁殖 并 随和 宿主 体液 流 问 全 
E, 同时 分 泌 能 抑制 或 影响 昆虫 免疫 与 代谢 功能 、 
扰乱 昆虫 赔 弃 和 变态 、 甚 至 致死 的 毒素 并 产生 草酸 
盐 结 唱 ， 导 致 生物 体 代 谢 雍 乱 ( Lewis et al., 2009; 
Ruiz-Sanchez ef al., 2010; Ruiz-Sanchez and O’ 
Donnell, 2012) ; 这 些 变 化 改变 了 家 看 对 一 些 具有 
特定 功能 重 日 的 需求 。 在 免疫 应 激 过 程 中 , AER 
被 重新 分 配 到 与 免疫 应 答 有 关 的 组 织 中 ， 以 合成 免 
疫 蛋 日 或 对 机 体 防 御 有 重要 作用 的 其 他 成 分 。 
此 , 受 刺 激 的 免疫 系统 扰乱 了 机 体 的 正常 生理 过 
程 , 也 诱导 了 其 对 一 些 特殊 氨基 酸 、 蛋 日 质 的 需求 
(Le Floc’h et al., 2004), Huang 等 (2009 ) 以 黑 胸 
败 血 病菌 感染 家 看 的 试验 也 证 实 部 分 基因 编码 的 重 
日 受到 调控 ,本 实验 中 抗原 人 处理 及 演示 通路 
( Trombetta and Mellman，2005 )、 阿 米 巴 病 通 路 
(Lejeune et al., 2009) 中 基因 表达 量 的 变化 与 该 说 
法 一 致 。 


核糖 体 通 路 中 的 核糖 体 蛋 白 与 机 体 细胞 的 增 
E FM. JE To HK Z8 2: UJ ( Lindström and Nistér, 
2010; Dutt et al., 2011) , 致使 其 差异 表达 与 机 体 的 
疾病 与 发 育 紧密 相关 , 机 体 病 变 时 核糖 体重 白 也 会 
低 水 平 表达 (Dai et al., 2007) 。 氮 酰 -tRNA 生物 合 
成 通路 中 的 氨 酰 -:RNA 连接 酶 既 能 专 一 性 地 识别 氨 
基 酸 的 侧 链 和 tRNA, 也 参与 了 转录 、 翻译 调 控 、 信 
号 转 导 、 凋 亡 及 其 他 免疫 应 激 等 过 程 (机 捷 和 金 由 
3E, 2004)。 凋 亡 也 是 昆虫 抗 感染 的 一 种 应 激 机 制 
(Feng et al., 2007) 。 剪 接 体 通路 中 , 剪接 体 的 剪接 
MIE mRNA 前 体 成 为 成 熟 、 有 功能 的 mRNA, 也 
是 形成 集 日 质 多 样 化 的 关键 机 制 之 一 ( Maniatis and 
Tasic, 2002) 。 本 研究 的 3 个 显著 性 富 集 的 通路 中 
相关 基因 表达 的 变化 说 明 , HERREZEN, Tu 
主 可 能 局 动 了 凋 亡 途径 , 合成 具有 免疫 功能 的 重 日 
来 抵御 其 体内 的 菌 丝 ; 另外 , 感染 后 宿主 的 食量 相 
对 减少 、 生 长 变 缓 , 导致 机 体 需要 合成 的 和 蛋白质 也 
相对 减少 , 致使 一 些 相关 基因 的 表达 水 平 下 调 。 

定量 PCR 分 析 表 明 , 输入 重 日 7 (importin 7, 
Imp-7 ) , 假定 的 溶菌 酶 (putative lysozyme) ~ 30kP E 
白 酶 A Hif (30 kP protease A precursor) 基因 的 表 
达 趋 势 与 数字 表达 谱 结 果 一 致 。Imp-7 属于 B 输入 
BAAR, 其 表达 水 平 的 变化 是 细胞 分 化 的 一 个 关 
PEJE ZS (Ao et al., 2007; Zaitseva et al., 2009) 。 
溶菌 酶 能 水 解 、 切 开 病 原 上 肽 聚 糖 的 B-1, 4-0 H AR, 
这 些 被 切 开 的 片段 由 于 能 被 昆虫 的 肽 聚 糖 识别 重 日 
识别 从 而 起 到 信号 分 子 的 作用 , 并 局 动 昆虫 的 天 然 
免疫 应 答 (Park et al., 2007) ; 溶菌 酶 还 可 调 广 外 源 
物 的 黑 化 (Li and Paskewitz, 2006), 以 细菌 为 食物 
来 源 并 清除 其 他 抗菌 蛋白 作用 于 细菌 时 产生 的 雄 片 
(Erban and Hubert, 2008) 。30kP 蛋白酶 A 是 类 胰 
肽 酶 丝氨酸 蛋白 酶 之 一 , 其 选择 性 水 解 的 30 kD Ez 
日 参与 了 酚 氧 化 酶 原 的 激活 并 能 抑制 看 体内 菌 丝 的 
生长 (Ujita et al., 2005). HF A tie ER ZZ FE A He 
AS. KARA I KENARA FATX 3 
HARD B FI PARE AA IB MIL, T6 
FS RS AS E ESD AY BUS NIT KERRIE 
化 或 消化 位 于 宿主 体内 的 菌 丝 ; 30kP EAR A 前 
体 可 能 促进 了 30kP 蛋白酶 A 的 转化 , 最终 促 使 特 
环 系统 中 30kD EA AY REM el ZER, M 
Imp-7 可 能 促进 了 相关 和 蛋 日 的 转运 。 

HK SS $F UH A (peptidoglycan recognition 
proteins, PGRPs ) 是 生物 体内 非常 保守 、 能 识别 病 
原 相 关 分 子 模式 的 一 类 模式 识别 受 体 ; 本 实验 中 
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PGRP-S6 是 一 种 低 分 子 量 的 短 链 胞 外 重 折 。PCRPs 
能 启动 昆虫 免疫 系统 中 酚 氧 化 酶 原 的 激活 ， 并 将 有 
毒 的 黑色 素 转移 到 入 侵 的 微生物 中 (Park et al., 
2006) 。 该 受 体 通常 作用 于 Toll 和 IMD 途径 的 上 游 
以 识别 、 连 接 病 原 体 壁 的 肽 聚 糖 (Tanji and Ip, 
2005), DGE 与 定量 PCR 方法 均 检 测 到 该 基因 的 表 
达 上 调 , 表明 PGRP-S6 识别 并 连接 到 真菌 的 肽 限 糖 
Es 启动 了 答 主 的 天 然 免 疫 。 

HERREZEN, 家 看 的 一 些 相关 基因 的 表 
达 发 生变 化 ; 其 他 病原 感染 家 看 后 ,也 有 相似 的 变 
化 出 现 。Zhou 等 (2013 ) HEN Fr 43 Pr 2€ AD. 
易 感 NPV 品种 感染 BmNPV 时 发 现 , 在 感染 0 ~ 48 
h 内 丝氨酸 重 白 酶 前 体 基 因 的 表达 量 在 所 有 家 看 品 
种 中 几乎 没 变 化 ; 而 在 感染 0 ~72 h 内 , 感性 品种 
中 的 表达 量 显 著 高 于 抗 性 品种 。 而 Qin 等 (2012 ) 
等 用 相同 的 试验 材料 得 出 : ZARE A M AR 
和 蛋白酶 前 体 的 转化 产物 ) 基因 在 抗 性 品种 中 显著 表 
达 ， 而 在 感性 品种 中 却 不 表达 。 目 前 依然 无 法 解释 
丝氨酸 重 日 酶 基因 的 表达 是 受 调 控 还 是 受 微生物 成 
分 的 诱导 。 尽 管 Sagisaka 等 (2010 ) 等 发 现 家 和 看 感染 
BmNPV 12 h 后 , ABE ARAN RADLER; 
Bao = ( 2009 ) 用 抑制 性 消减 杂交 技术 、Wu 等 
(2011) 等 用 芯片 技术 却 发 现 家 不 感染 BmDNV 及 
BmCPV 5, ABE A ZEA ZA Eva; 本 实验 中 也 检 
Tui BUY AS A (Hspl, Hsp21. 4, Hsp25. 4 和 Hsp60 
等 ) 基 因 的 表达 上 调 。 这 些 说 明 不 同类 型 的 外 源 物 
在 不 同 的 家 看 品种 间 会 导致 应 党 反应 的 差异 。 

本 人 研究 设计 了 10 个 差异 表达 基因 的 灾 光 定量 
PCR 引物 进行 定量 PCR 实验 以 验证 数字 表达 谱 的 
准确 性 ; 结果 发 现 其 表达 趋势 一 致 (图 5$) 。 由 于 用 
不 同 的 研究 方法 得 到 的 基因 表达 结果 可 能 不 一 致 
(Luo et al., 2005) ,所 以 用 数字 表达 谱 获 得 的 差异 
基因 表达 状况 需要 用 定量 PCR 验证 。 

本 研究 用 高 通 量 的 Solexa 测序 技术 对 感染 白 从 
菌 的 家 和 看 幼虫 进行 转录 组 测序 和 功能 分 析 , RAR 
得 5 472 个 基因 , 其 中 有 377 个 差异 表达 基因 ，12 
个 上 调 基 因 、26 个 下 调 基因 参与 3 个 显著 性 宦 集 的 
KEGG Pathway。 但 关于 这 些 基 因 的 差异 表达 及 其 
参与 的 Pathway ERAH H E A IS zd Te 9 RE 
与 作用 , 需要 今后 进一步 深入 开展 人 研究。 家 看 感染 
日 仅 菌 的 转录 组 研究 , WARS ABA ELTE AR 
看 各 种 途径 相关 的 关键 基因 的 发 掘 为 基因 克隆 、 功 
能 研究 提供 了 数据 基础 ， 为 家 和 看 抵御 真菌 侵 染 的 研 
究 、 日 僵 阔 应 用 于 农林 害虫 的 生物 防治 新 策略 的 人 研 


发 提供 了 丰富 的 基因 信息 及 新 的 研究 方 癌 。 
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